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Abstract 

We have performed numerical simulations of the kinematic induction equation in order to examine the dynamo efficiency 
of an axisymmetric von-Karman-like flow subject to time-dependent non-axisymmetric velocity perturbations. The numerical 
■ model is based on the setup of the French Von-Karman-Sodium dynamo (VKS) and on the flow measurements from a model 
' water experiment conducted at the University of Navarra in Pamplona, Spain. 

Our simulations show that the interactions of azimuthally drifting flow perturbations with the fundamental drift of the 
magnetic eigenmode (caused by the inevitable equatorial symmetry breaking of the basic flow) essentially determine the 
temporal behavior of the dynamo state. We find two distinct regimes of dynamo action that depend on the (prescribed) drift 
frequency of an (m = 2) vortex-like flow perturbation. For comparatively slowly drifting vortices we observe a narrow window 
with enhanced growth-rates and a drift of the magnetic eigenmode that is synchronized with the perturbation drift. The 
j— | ■ resonance-like enhancement of the growth-rates takes place when the vortex drift frequency roughly equals the drift frequency 
of the magnetic eigenmode in the unperturbed system. Outside of this small window, the field generation is hampered compared 
to the unperturbed case, and the field amplitude of the magnetic eigenmode is modulated with approximately twice the vortex 
drift frequency. The abrupt transition between the resonant regime and the modulated regime is identified as an spectral 
exceptional point where eigenvalues (growth-rates and frequencies) and eigenfunctions of two previously independent modes 
collapse. 

i-C . In the actual configuration the vortex drift frequency observed in the water experiment is much larger than the fundamental 

drift frequency of the magnetic eigenmode obtained from simulations and we conclude that the fulfillment of the resonance 
condition might be unlikely in present day dynamo experiments. Nevertheless, these radial vortices can be anchored to fixed 
I . spatial azimuthal positions (using holes or fingers) breaking the axisymmetry of the experimental setup. Moreover, a parametric 
^ • resonance might occur in the Earth's core where tidal forcing provides appropriate flow perturbations, or in future dynamo 
y—{ ' experiments that will allow a better control of periodic non-axisymmetric velocity modes in order to meet the resonance 
\f} , condition. 



m 

O 

(N 



- 1—1 

x 



PACS numbers: 47.65.-d, 91.25.Cw, 52.65.Kj 

Keywords: kinematic dynamo simulations, dynamo experiment, parametric resonance 



* a.giesecke@hzdr.de 



1 



I. INTRODUCTION 



Cosmic magnetic fields are ubiquitous phenomena that are intrinsically coupled to most astrophysical objects like 
planets, stars or galaxies. The origin of these fields involves the formation of electrical currents by means of a complex 
flow of a conducting fluid or plasma. This process, the so called dynamo effect, is necessarily three dimensional 
and non-linear which makes an analytical or numerical approach difficult. Meanwhile, fluid flow driven generation 
of magnetic fields has also been obtained in laboratory experiments providing a complementary tool to astronomical 
observations or numerical simulations. However, whereas astrophysical dynamo action is comparably easy because of 
the large dimensions of the involved flows, its experimental realization requires considerable technical efforts An 
important obstacle for the occurrence of laboratory dynamo action arises from the scaling behavior of the power that 
is required to drive a flow with a requested magnetic Reynolds number Rm. For turbulent flows this power scales 
cx Rm 3 so that a reduction of the critical Rm for the onset of dynamo action is most important to achieve magnetic 
self-excitation at all. 

So far dynamo experiments based on a flow of a conducting fluid have been successfully conducted in Riga Q, 
Karlsruhe [|[ and Cadarache Q . The first two facilities made use of a more or less predetermined fluid flow essentially 
fixed by the forcing and the shape of the internal tubes. Note, however, that at least in the Riga dynamo experiment 
the saturation process involved already a non-trivial back-reaction effect of the magnetic field that changes the 
geometry of the flow. Such effects might be even more pronounced in the Cadarache von-Karman-sodium (VKS) 
dynamo. In that experiment the flow driving by two opposing impellers provides more freedom for the development 
of a saturated turbulent state, in which the back-reaction of the magnetic field on the fluid can strongly modify the 
geometry and dynamics of the flow. In an idealizing model the mean axisymmetric flow between counter-rotating 
impellers comprises two toroidal and two poloidal eddies (so called s2t2 topology) , and it is well known that this flow 
is able to drive a dynamo @, Q . Various attempts in different geometries have been made (numerically as well as 
experimentally) in order to examine dynamo action driven by such a flow [4, 7-20]. However, so far experimental 
dynamo action driven by a von-Karman-like flow is obtained only at the VKS facility and only when at least one of the 
flow driving impellers is made of soft iron with a large relative permeability. Kinematic simulations of the Cadarache 
dynamo demonstrated the close linkage between the exclusive occurrence of dynamo action in the presence of soft iron 
impellers and the observed axisymmetry of the magnetic field [2lll22l|. Nevertheless, a fully satisfactory explanation 
of the working principle of this dynamo is still missing and it is still unclear whether the present experiments will ever 
be able to achieve growing equatorial dipole modes, which constitute the originally expected magnetic field geometry. 

An improvement of present numerical models may require the explicit consideration of coherent non-axisymmetric 
structures that repeatedly have been observed in water experiments using a von-Karman-like flow driving [23T4251 ] . 
Non-axisymmetric time-periodic flows with a dominant azimuthal wavenumber m — 2 have also been found in 3D 
simulations of s2t2 flows in spherical geometry [lj| . Kinematic dynamo simulations using various manifestations of 
these velocity fields showed a surprising diversity of behavior patterns: self-generation of magnetic energy was only 
found when the time-dependent flow field was taken into account whereas the simulations with the time averaged flow 
or with different snapshots of the velocity field did not exhibit dynamo action. 

The present paper aims at a numerical examination of the dynamo efficiency of a cylindrical VKS-like system 
subject to non-axisymmetric velocity perturbations with a single azimuthal wavenumber m = 2. Such velocity modes 
were observed in a water experiment conducted at University of Navarra in Pamplona in order to analyze the influence 
of slowly evolving large scale flow on the occurrence of dynamo action [23], (2(| . In this study we utilize the essential 
features of the measured flow field as the basis input for numerical simulations of the kinematic induction equation. 
Typical input parameters that are systematically varied are the flow amplitude (in terms of the magnetic Reynolds 
number) and the azimuthal drift motion of the implied non-axisymmetric velocity perturbation. From the simulation 
data we extract the leading eigenmodes and the related eigenvalues in terms of growth-rates and frequencies that 
describe field amplitude modulations and/or azimuthal field drift. Interestingly, for comparably low drift frequencies 
of the velocity perturbation we first observe a phase locking of the magnetic eigenmode drift with the vortex drift 
which is replaced, for higher perturbation drift frequency, by the appearance of a time modulated magnetic eigenmode. 
By analyzing the involved growth rates and frequencies in the phase-locked regime and in the modulated regime, we 
identify the transition between them as a spectral exceptional point where eigenvalues and eigenfunctions of two 
modes coincide [27|, [28|. The observed behavior is in close analogy with typical (resonant) mechanical systems subject 
to periodic forcing, like e.g. spinning disk systems (2i| or to the behavior observed in the stability study of water 
waves 30] . Comparable effects have also been found in mean-field dynamos of aw-type that were designed to explain 
the bisymmetric field pattern observed in spiral galaxies. In these models a periodic perturbation is caused by density 
waves due to spiraling arms, and a parametric resonance (also called swing excitation) is observed when the frequency 
of the perturbing velocity pattern is twice the oscillation frequency of the (axisymmetric) dynamo (3l| - [35| . 
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II. EXPERIMENTAL BACKGROUND 



The simulations presented below are motivated from a water experiment that is described in detail in [23[. In 
the experiment a von-Karman swirling flow is driven by two counter-rotating impellers located at the end caps of a 
cylindrical vessel of radius R — 10 cm and height H = 20 cm ( figure QJi). Both impellers spin with a rotation rate up to 




(a) (b) (c) 

FIG. 1. (a) Sketch of the experimental setup. In the experiment the von-Karman swirling flow is driven by two counter-rotating 
impellers located at the end caps of a cylindrical vessel of radius R — 10 cm and height H — 20 cm. The yellow ring denotes the 
equatorial plane of the cylinder, (b) Instantaneous snapshot of the turbulent fluid flow. Note the vortex like structure roughly 
in the center of the vessel, (c) Mean axisymmetric flow field measured in the experiment. The colored structure denotes the 
toroidal flow (u v ) and the arrows denote the poloidal flow (u r ,u z ). Although both impellers rotate with the same velocity 
different sizes of the flow cells occur due to spontaneous symmetry breaking. 

12 Hz so that the resulting flow is highly turbulent with flow fluctuations of the same order as the mean flow (a typical 
snapshot of the turbulent flow is shown in figure [1}d). The velocity field is measured using Laser Doppler Velocimetry 
(LDV) and the mean velocity field is obtained by averaging the instantaneous flow for (at least) 100 impeller turns. 
The resulting axisymmetric velocity field consists of two toroidal cells and two meridional recirculating cells that 
are roughly restricted to each cylindrical half-space (figure [TJ;) . Comparable structures have been also observed in 
water experiments with similar configurations [24l |25| . The cylindrical configuration combined with the specific flow 




(a) (b) 

FIG. 2. Reconstruction of averaged trajectories in the experiment. The red and blue lines represent the averaged trajectories 
impelled by each one of the top and bottom impellers, (a) Symmetric flow configuration, (b) Equatorial symmetry breaking 
observed in the experiment when Re > 10 4 . Spontaneous jumps between both possible states are observed in the experiment. 

driving imposes a discrete symmetry: the most symmetric flow that can be obtained between two propellers rotating 
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in opposite directions breaks the reflection symmetry on any plane that contains the cylinder axis, but preserves 
the symmetry about a 7r-rotation around any diameter in the equatorial plane (figure [2^). However, in the water 
experiment this symmetry is broken for Re <; 10 4 (without internal symmetrizing fixtures) and the measured mean 
velocity field becomes asymmetric even when both impellers spin with equal rotation rates ^23j. As a consequence 
one of the cells (red or blue in figure [2}d) becomes larger than the other. The averaged velocity field is meta-stable, 
and spontaneous jumps of the dominant cell from one side to the other are observed (so called inversions) with both 
possible states occuring with equal probability. 

The transition time for an inversion is Ti nv ~ 10 s and the time between consecutive events varies from minutes to 
hours following a probability distribution p(t) = 1 /T e~ t / T ° with a parameter To determined by the noise intensity. Tq 
decreases for increasing Reynolds number, i.e. inversions occur more frequently if the turbulence is more vigorous. 
The inversions of the velocity field introduce a slow time scale at which the large scale Sow and hence the large scale 
magnetic field change their structure. In the long run, the enhanced dissipation arising from the transition between 
different eigenstates (corresponding to the different velocity fields) can result in a reduced efficiency for the dynamo 
process when the mean life time of a velocity state becomes of the same order as the magnetic diffusion time [361 ] . 

In addition to breaking of the equatorial symmetry the observed mean flow also violates the ideal axial symmetry. 
Coherent non-axisymmetric structures emerge close to the equatorial plane (even before a turbulent state is reached) 
and establish a local swirling flow around an axis perpendicular to the main symmetry axis of the cylinder. The non- 
axisymmetric flow perturbation is dominated by an azimuthal wavenumber m = 2 and undergoes an azimuthal drift 
that is immediately linked to the flow orientation of the dominating toroidal cell. The drift frequency of the vortex- 
like pattern is related to the maximum azimuthal mean flow velocity it™ ax by w v ~ 0.3u™ ax /i? v , where R v denotes 
the radius of the maximum vortex velocity (R v w 0.857 in units of the cylinder radius). The mid-size structures 
provide an additional source of helicity, so that it is likely that the vortices are involved in the dynamo process when 
a conducting fluid like liquid sodium is utilized. 



III. NUMERICAL MODEL 

The temporal development of the magnetic flux density B induced from a flow u of a conducting fluid is described 
by the magnetic induction equation that results from the combination of Faraday's law, Ohm's law and Ampere's law 
(without displacement current): 

— = V x (u x B - nV x B) . (1) 

In ([T]) 7/ denotes the magnetic diffusivity which is related to the electrical conductivity a and the vacuum permeability 
Ho by n = (po°~)~ 1 - Equation ([T]) is time stepped applying a finite volume method where a constraint transport 
scheme ensures the exact treatment of the solenoidal property of B (if the initial field is divergence free) . Insulating 
boundary conditions are treated with a modified boundary integral equation approach which yields the tangential field 
components on the boundary from the normal field components on the whole surface of the computational domain 

[Ull. 



A. The axisymmetric velocity field 

We use an analytically prescribed flow field that incorporates the main characteristics of a von-Karman flow and 
allows a convenient variation of the equatorial symmetry breaking. Since the typical life-time of the mean flow state 
is much longer than the time scales governing the (non-axisymmetric) flow perturbations the dynamics introduced by 
the inversions can be ignored so that we assume a time-independent axisymmetric velocity field. Equatorial symmetry 
breaking is modeled using a basic flow that is composed of two different axisymmetric parts with different symmetry 
properties. The main part is symmetric with respect to a 7r-rotation around any diameter in the equatorial plane. 
This contribution will be called even flow (figure 13^.). The geometric structure of the even flow consists of two counter- 
oriented toroidal cells and two recirculating poloidal cells located in each cylindrical half-space. Mathematically, the 
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FIG. 3. Prescribed velocity field. From left to right: (a) even flow (MND flow) as given by equation (J2j) , (b) odd flow 
contribution given by equation ([3]), (c) total flow u — u c + au° for a — 0.35. The color coded pattern represents the azimuthal 
flow {u v ) and the arrows denote orientation and magnitude of the meridional flow (u r ,u z ). 



even part is prescribed by the so called MND-flow first proposed by Marie, Normand and Daviaud |39j |: 



u°(r,z) = -— r(l -r) 2 (l + 2r)cos 



H 



u e v (r, z) = 4er(l — r) sin 



(2) 



u° z (r,z) = (1 — r)(l + r — 5r 2 ) sin 



2ttz 



where e represents the relation between toroidal and poloidal flow (here e = 0.8155) and H is the total height of the 
cylinder (here H = 2.0). 

The second contribution changes its sign when performing the same 7r-transformation and is called odd flow (fig- 
ure [3Jd). This contribution is determined by a single poloidal cell that spreads through the whole cylinder and a global 
azimuthal rotation u v (independent of z and ip). The odd flow contribution assumes the same expressions for e, H 
and for the radial profile as used for the even flow and is given by: 



-r(l - r) 2 (l + 2r) sin 



2lTZ 



' 2H 

u° = (l-r)(l+r- 5r 2 )cos 2 (J^j . 



(3) 



The variation of the equatorial symmetry breaking is realized by multiplying the odd flow contribution ([3]) with a 
weighting factor a £ [0; 1] so that the total flow field is given by 



u = u 



au 



(4) 



u describes a flow field with two different cell sizes with the even contribution always being dominant (see figure |3j;) . 
In our model the dominant toroidal eddy is always located close to the north (top) impeller (at z = +H/2) and is 
characterized by a counterclockwise oriented azimuthal velocity (figure [5t). 

The combined flow field is scaled allowing a systematically variation of the (prescribed) magnetic Reynolds number 
defined as 



Rm 



WZ 

V 



(5) 



where rj denotes the magnetic diffusivity, 1Z a characteristic length scale (here 1Z = 1, the radius of the cylindrical 
domain) and U is the peak velocity defined by U — U max = max[(u 2 + w 2 + m 2 ) 1 ^ 2 ]- 

The flow field (U) is qualitatively similar to the model flow applied in (4p| where dynamo action is examined in a 
VKS-like configuration with impellers rotating at different speeds. There are various other possibilities to realize the 
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equatorial symmetry breaking of the ideal flow field ([2]) and it turns out that not every property of the observed flow 
field can always be reproduced exactly. For example, the definition of the flow field (|4|) does allow an easy adjustment 
of the symmetry breaking in terms of the parameter a, but exhibits different azimuthal velocities near the impellers. 
In this sense the idealized flow Q differs from the observed flow which results from a forcing through impellers that 
are in exact counter-rotation, so that the flow velocities close to both impellers are the same. Such deviations slightly 
influence the quantitative outcome of the simulations like e.g. the critical magnetic Reynolds number for the onset of 
dynamo action, nevertheless they do not affect the essential conclusions that will result from our simulations. 



B. Modeling of the non-axisymmetric velocity perturbations 

The coherent non-axisymmetric flow perturbations are located near the outer wall of the cylinder and close to 
the equatorial plane where a strong shear layer emerges that is caused by the opposite azimuthal flow orientation in 
each cylindrical half. The vortices allow a relaxation of the shear in the equatorial layer in some way similar to a 
Kelvin-Helmholtz instability, but under turbulent conditions. The formation of the vortices in the experiment and 
the corresponding implementation in the numerical model are sketched in figure [U In the simulations the vortices 




FIG. 4. Mean flow and vortex pattern realized in the model: The figures show axisymmetric and non-axisymmetric contributions 
of the velocity field in the (tp, z) plane at r — 0.9. Only the azimuthal and axial components of the velocity field are presented, 
(a) odd and even axisymmetric part (eqs. ([2]) and (|3|). (b) Vortices according to (|6} for m = 2. (c) Total velocity field 
(combining eqs. ([2]), (fSJl and <j6j ) for V m = 0.2. The colored structures represent the radial component of the vorticity of the 
non-axisymmetric contribution (V x u v ) r (red (blue) corresponds to positive (negative) vorticity). (d) Same as in (c) for the 
experimental data. 

are parameterized by analytical expressions whereby the essential properties - the azimuthal wave number, thickness, 
diameter and average radial position are taken from the measurements. The center of the vortices with a typical 
diameter of 0.5 (and an aspect ratio of unity) is located around the equatorial plane with a maximum at a radius of 
R v = 0.857 (all numbers are denoted in units of the cylinder radius R = 10 cm). The non-axisymmetric contribution 
of the vortices to the velocity field is explicitly given by 

u v r =0 

ul = — cos(m(<p + cj v t))r 2 (l - r )e-^ r - ro ^' sm (4 z ) (6) 
u v z = -2V m sm(m(ip + io v t))r 2 {l - r )e-^ r - ra ^^ cos 2 (2z), 

with the axial coordinate z restricted to the interval z £ [— 7r/4; 7r/4]. In © Vm denotes the magnitude of the 
non-axisymmetric perturbation and the quantities a and ro are estimated from the experimental observations (here 
always ro = 0.9 and a = 0.12). The azimuthal drift frequency of the vortex perturbation is denoted by u> v , which 
in the experiment is linked to the azimuthal velocity of the dominant mean flow cell. In the system of units applied 
in the simulations the observed vortex drift frequency ~ 0.3u™ ax /i? v corresponds to ui v w 17.1 (14.9, 14.0, 13.1) for 
a = 0(0.62,0.84,1.00). Note that the drift frequency of the non-axisymmetric velocity pattern with an azimuthal 
wave number m is related to the previously defined vortex frequency by 57 p = mw v . In the water experiment the 
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vortices exhibit an intermittent behavior, e.g. sudden jumps to a state with a different wavenumber (m = 3 or m — 4) 
or fluctuations of the vortex drift. Furthermore, the locations of the vortices undergo variations that are connected 
to the largest fluctuations of the mean velocity field. Those fluctuations appear centered around the shear layers, 
and close to the wall. These temporal alterations are not considered in our numerical model, where we restrict our 
examinations to the case m = 2 which corresponds to the most probable configuration observed in the experiment. 
In all cases the disturbance introduced by the non-axisymmetric component remains only weak and has no influence 
on the actual Rm given by the maximum of the modulus of it. 



IV. RESULTS 

In the kinematic approach the back-reaction of the magnetic field on the flow is neglected so that the solution of 
the induction equation ([1]) represents a linear problem, which in principle could be solved with the Ansatz B(r,t) = 
B (r)e xt . A dynamo solution is obtained if the magnetic energy density E mag = (2^ ) _1 |-B| 2 grows exponentially 
oc e 2A *. In general A is a complex quantity A = 7 + iu where 7 denotes the field amplitude growth rate and to denotes 
an oscillation or drift frequency In accordance with Cowling's anti-dynamo theorem, dynamo solutions generated by 
a prescribed mean axisymmetric flow necessarily yield a non-axisymmetric field. In all simulations presented below 
the (m — 1) mode is the dominant field contribution so that the magnetic eigenmode behaves oc cos<^ and represents 
a dipole oriented along the equatorial plane. This is in contrast to the real VKS experiment in which the presence 
of the magnetic impellers works in favor of an axisymmetric eigenmode (axial dipole, [2ll. |22| ). In the following the 
growth rates represent magnetic field amplitude growth rates for this (m — 1) mode and the applied time scale is 
given by the magnetic diffusion time t v — R 2 /rj. 



A. Equatorial symmetry breaking 

For a sufficiently large magnetic Reynolds number the flow field Q drives a dynamo with a typical structure shown 
in figure [5] which displays the eigenmodes for the undisturbed flow field (a — 0) and the model case a = 0.62. In 
both cases the geometry is dominated by two interleaved banana-cell like structures. However, in case of equatorial 
symmetry breaking (a = 0.62, right panel) a slight concentration of magnetic energy is observed in the upper cylindrical 
half-space containing the dominant flow cell. The breaking of the ideal equatorial symmetry suppresses dynamo action: 
when a symmetry breaking flow contribution is added the critical magnetic Reynolds number for the onset of dynamo 
action roughly increases oc a 2 from Rm c ss 50.5 for a — to Rm c = 109.5 at a = 1 (figure [6^ & b). 

In the ideal symmetric case (a = 0) the eigenmode shows no time dependence except the exponential decay. 
With symmetry breaking the (m = 1) mode exhibits an azimuthal drift around the cylinder axis, i.e. B ~ 



FIG. 5. Geometric structure of the dynamo eigenmode. Left: no equatorial symmetry breaking (a — 0), Rm = 70. Right: 
a — 0.62, Rm = 70. The isosurfaces show the magnetic energy density at 10%, 30% and 50% of the respective maximum value. 
The colored contours denote the azimuthal magnetic field B v . Note the magnetic energy concentration in the upper part of 
the cylinder which contains the dominant flow cell in case of equatorial symmetry breaking (right panel). 
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FIG. 6. (a) Growth-rate versus Rm for various equatorial symmetry breaking a; (b) critical magnetic Reynolds number versus 
a. The dashed curve denotes a fit by a polynomial of degree 2. 



Bo(r, z)e xt cos(ip — LUft) with the magnetic field drift frequency Wf increasing nearly linearly with the degree of asym- 
metry a and with Rm (figure [7^, & b) . u>{ is constant in time and is oriented opposite to the azimuthal flow in the 
dominant cell. For axisymmetric velocity fields (as u e and u°, eqs^l&E]) poloidal and toroidal velocity components 
are decoupled, and their shear layer may have independent locations. The relative position of these shear layers 
controls the magnitude (and orientation) of the magnetic field drift. For example, an exclusive increment of the odd 
azimuthal flow contribution (which only shifts the toroidal shear layer) would alter the drift frequency. The change 
might be quite dramatic and even a reversal of the drift orientation is possible with the eigenmode co-rotating with 
the dominant azimuthal flow, as observed for example in [8j. A systematic study of the influence of this effect is under 
way. 



(o) (b) 




FIG. 7. (a) Drift frequency of the (m = 1) mode versus Rm for various a. (b) drift frequency of the (m = 1) mode versus 
equatorial symmetry breaking a. Note the slight deviation from the linear behavior for small a. 



B. Impact of non-axisymmetric velocity perturbations 

The presence of non-axisymmetric velocity components results in a coupling between different azimuthal magnetic 
field modes that are separated in the purely axisymmetric problem. Hence the characterization of the dynamo 
eigenmode with a single azimuthal wavenumber is no longer correct because (in principle) all azimuthal wave numbers 
are linked. In the following, we examine the impact of a non-axisymmetric flow perturbation with a (single) wave 
number m = 2 so that two classes of magnetic eigenmodes arise which incorporate even azimuthal wave numbers 
(to = 0,2,4,...), or odd azimuthal wave numbers (m = 1,3,5,...), respectively. In our model only the second class 
with odd wave numbers is relevant whereas the "even modes" decay on a much faster timescale. 
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In the following, we keep the amplitude of the non-axisymmetric velocity contribution fixed at V m = 0.3C/ max which 
roughly corresponds to the value observed in the experiments. In this regime the presence of the non-axisymmetric 
velocity component neither changes the global (axisymmetric) flow topology nor the actual magnetic Reynolds number. 
We further assume that the drift frequency of the vortex, cj v , is a free parameter that is systematically varied in the 
interval uj v £ [—60, +60]. The maximum/minimum values correspond to a vortex drift velocity approximately equal 
to the maximum azimuthal flow velocity (recall that in the water experiment the observed vortex drift velocity is 
roughly 0.3u™ ax ). We also examined the artificial case with a vortex drift orientation opposite to the azimuthal flow 
of the dominant cell. These cases are denoted by negative frequencies. The temporal development of the amplitude of 




E 
< 




FIG. 8. (a) Time dependence of the magnetic field amplitude for Rm = 70 and a = 0.62. (b) Fourier spectrum for the case 
cj v = 14.9 (after eliminating the exponential decay). The various peaks denote the modulation frequency of the field amplitude 
w a re 25.08 (and its overtones). 

the dominant eigenmode for two different vortex drift frequencies, but similar Rm, is shown in figure [5] (dashed curve: 
w v re 2.5, solid curve: cj v re 14.9). Both curves show completely different behavior: for lo v re 14.9 the solution decays 
and its amplitude is modulated with a frequency of w a re 25.08 (solid curve), whereas at ui v re 2.5 the field amplitude 
does not show any time dependence except the exponential growth (dashed curve). Regardless of the distinct temporal 
behavior for both cases the field structure is quite similar (except the pulsating character in the modulated case) and 
remains close to the pattern already observed in the previously discussed unperturbed configuration (figure [9]). 





FIG. 9. Left: field structure in the resonant regime (Rm = 70, a = 0.62, u v = 2.5, corresponding to the dashed curve in the 
left panel of figure [SJ). Right: time series of a typical solution with amplitude modulation (Rm = 70, a — 0.62, u v re 14.9, 
corresponding to the solid curve in the left panel of figure[8]|. The isosurfaces show the magnetic energy density at 10%, 30% and 
50% of the respective maximum value and the colored contours denote the azimuthal magnetic field B v . The time series covers 
one modulation period with each snapshot scaled by the respective maximum so that the effects of exponential decay/growth 
and the amplitude modulations are eliminated. The period of the amplitude modulation is much shorter than the time scale 
of the magnetic field drift so that the phase of the eigenmode (i.e. its orientation in space) remains nearly constant. 
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FIG. 10. (a) Growth-rate versus vortex drift frequency for Rm = 70 and V m = 0.3C/ max . In the abscissa the resonance frequency 
CJv SS has been subtracted so that all curves are centered at the origin. The horizontal dotted lines denote the growth rates 
without non-axisymmetric perturbation, (b) enlarged section of the resonant regime. The growth-rates in the resonant regime 
correspond to non-oscillating m = 1 eigenmodes whereas outside of this window the solutions exhibit a modulation of the field 
amplitude (see also figure \TTk ). 



The occurrence of non-axisymmetric solutions with oscillating energy is a clear indication for the presence of two 
distinct azimuthal modes with the same growth-rate (i.e. real part of the eigenvalues) but with different frequencies 
(imaginary part of eigenvalues). In fact simulations that cover a sufficient number of periods exhibit characteristic 
beat patterns indicating that the frequencies of the superimposed modes are very close. We observe field amplitude 
modulations independently of the degree of equatorial symmetry breaking for a broad range of vortex drift frequencies, 
whereas unmodulated solutions only occur within a narrow window of relatively slow vortex drift frequencies. The 
unmodulated solutions are further characterized by a sharp maximum for the corresponding growth-rates (figure HOD , 
typical for parametric resonance. Characteristic properties of the resonant regime and a comparison to values in 
the unperturbed state are denoted in table |TJ The location of the resonance maximum (denoted by lu^ cs ) is roughly 



TABLE I. Characteristics properties of the resonant window for different values of the equatorial symmetry breaking a. All 
data stems from runs with Rm = 70. The critical magnetic Reynolds number (column 6 & 7) has been estimated from 
inter/extrapolation of simulation runs with Rm = 50, 70 and 90. 
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growth rate 
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unperturbed 
growth rate 




Rm c 
without 
vortex 


Rm c at 
reso- 
nance 


relative 
reduction 
of Rm c 


width of 
unmodulated 
regime Aoj v 


0.00 


+2.817 


-0.00 


+0.912 


0.00 


59.7 


50.5 


15.4% 


5.06 


0.62 


+0.667 


-1.99 


-0.798 


-2.63 


78.5 


64.8 


17.5% 


5.11 


0.84 


-0.302 


-2.91 


-1.839 


-3.41 


92.1 


72.9 


20.8% 


5.63 


1.00 


-1.181 


-3.50 


-2.905 


-4.04 


108.5 


82.7 


23.8% 


6.12 



determined by the field drift frequency of the magnetic eigenmode in the unperturbed case (from now labeled with ui®). 
However, a constant gap is observed in runs with equatorial symmetry breaking (a > 0) resulting in lj™ s ~ (wf +0.5), 
whereas the case a = remains a particular case with w™ s = 0. The shift of the resonance frequency with respect 
to the fundamental frequency in the unperturbed problem probably follows from the interaction of vortex structure 
and drifting eigenmode similar to the shift of the resonance frequency in periodically forced mechanical systems with 
damping. Outside of the resonance, the growth-rate is suppressed compared to the unperturbed case (denoted by the 
dotted horizontal line in figure [T0| . Only for very large drift velocities the value of the purely axisymmetric case is 
obtained again, i.e. the impact of the vortices vanishes when their drift frequency becomes too large. The width of 
the resonant regime slightly increases with the symmetry breaking parameter a and so does the relative enhancement 
of the growth-rates and accordingly the relative reduction of the critical Rm. 

The periodic distortion caused by the drifting vortices also influences the drift behavior of the magnetic eigenmode. 
In the resonant regime the dynamo eigenmode exhibits an azimuthal field drift that is immediately locked to the vortex 
drift frequency, i.e. the magnetic field pattern follows the vortices (wf ~ ui v , figure ITTb ) . This frequency locking also 
applies without equatorial symmetry breaking (a = 0) which does not show any field drift without vortices (red curve 
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FIG. 11. (a) average field drift frequency versus the vortex drift frequency for various values of the symmetry breaking parameter 
a. The solid curves denote the average field drift in the modulated regime. The horizontal dotted lines denote the field drift in 
the unperturbed case (wp ). (b) field drift frequency in the resonant regime. The field drift is synchronized with the vortex drift: 
tjf ~ Uy. (c) field amplitude modulation for various values of the symmetry breaking parameter a. The amplitude modulations 
scales cj a ~ 2w v . 



in figure [Tib). Beside the temporal decay/growth, the (constant) azimuthal field drift determines the only time scale 
in the resonant regime, whereas outside of the resonant window two different time scales appear: an average azimuthal 
field drift that is roughly determined by the equatorial symmetry breaking parameter a (solid curves in figure Hia tal 



This does not hold in the vicinity of the transition between the modulated regime and the resonant regime where the field drift exhibits 
a divergent trend. 
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FIG. 12. (a) compilation of frequencies of the magnetic eigenmode versus the vortex drift frequency minus the resonance 
frequency. The solid curves show a combination of field drift and field amplitude modulation and the dashed central part 
denotes the field drift locked to the vortices in the resonant regime, (b) corresponding growth-rates, Note the dashed curve 
in the bottom part which represents the assumed growth-rates of a hidden mode that cannot be identified with the applied 
time-stepping approach. 



and the aforementioned modulation of the field amplitude which scales ~ 2cj v (figure ITTb). In the modulated regime 
the field drift is not constant but varies regularly with the same timescale as the amplitude modulation. A movie that 
shows the time development of four characteristic runs (no vortex, resonant run with frequency locking and modulated 
case with positive and/or negative vortex drift frequency) can be found in the supplemental material for this article 

m- 

In order to illustrate the behavior of the eigenmode in the vicinity of the transition between the resonant and the 
modulated regime we have combined the curves shown in figure [11] into one drawing. Figure Q2] shows a combination 
of frequencies of the magnetic eigenmode (cjf + w a and ujf — oj a , left panel) and the related growth-rates (right panel) 
against uj v — w™ s . By using the sum/ difference of the observed frequencies (drift and amplitude modulation) a 
continuous transition but with a jump in the derivative is achieved between modulated and resonant regime. For 
larger vortex drift, the behavior is dominated by the amplitude modulation, whereas close to the transition (where 
w a — > 0) the behavior is determined by the field drift frequency ujf . 

The abrupt transition from the modulated regime with w a ~ 2uj v to an unmodulated regime with u a = and ujf ~ uj v 



occurs when two different eigenmodes merge. At these points, known as exceptional points [2jj, the eigenvalues of two 
eigenmodes coincide. In our case the change of the temporal behavior of the eigenmode results from the merging of 
the imaginary parts of the two eigenvalues whereas the real parts (growth-rates) are identical in the modulated regime 
giving rise to the field amplitude modulation. In the resonant regime (where the field amplitude modulation vanishes) 
the real parts of the eigenvalues of both interacting modes split and presumably form a "bubble" similar to the 
behavior of periodical perturbed resonant mechanical systems (see discussion below). In addition, at the exceptional 
point the two previously linearly independent eigenfunctions collapse and become indistinguishable. Mathematically 
this is described by the formation of a non-diagonal Jordan block structure in the matrix representation of the (non- 
selfadjoint) dynamo operator associated with algebraic eigenvectors [42j ■ A characteristic property of exceptional 
points with two coinciding eigenvalues is reflected in the time dependence of the field amplitude which, exactly at the 
exceptional point, exhibits an additional secular term linear in t [43{ so that 

B(t)^(b 1+ b 2 t)e^. (7) 



Such a time dependence has been found experimentally in an examination of a two level system in microwave billiards 
[lij and it would be interesting to identify this term in our simulations. However, up to present, the determination 
of growth-rates and frequencies near the transition from the modulated solution to the resonant solution is difficult 
and did not allow a reliable and unique determination of such a state. Hence, we can only speculate on the nature of 
the interacting eigenmodes because, unfortunately, with our present time-stepping method, we can only identify the 
leading eigenmode, whereas the identification of the second mode would require the use of an appropriate eigenvalue 
solver. However, it is rather suggestive to assume that the leading mode essentially consists of a dominant (m = 1) 
component with some slight addition of a (m — 2)-vortex induced (m = 3) component, and that the second eigenmode 
is dominated by a (m = 3) component, to which a slight (m = 2)-vortex induced (m = 1) component is added. 
Independently, we expect that the growth-rates of this "hidden" mode qualitatively behaves in the same way as in 
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comparable mechanical systems forming a bubble in the resonant regime. The assumed development has been added 
in figure [12] in terms of the dashed orange curve. This behavior is known e^g. from mechanical systems subject to 
periodic perturbations, the stability analysis of water waves conducted in [30] or the response of a non-linear state 
to spatially periodic forcing examined in (45j . A comparable pattern is also observed in the eigenvalue spectrum of 
oscillating a 2 -dynamos with non-trivial radial distribution of a 0,113]. However, in those models the behavior of the 
growth-rates and the frequencies is exchanged: the real part of the eigenvalues merge/split (in the same manner as the 
frequencies in the left panel of figure I12[) and the imaginary parts form a bubble corresponding to a restricted regime 
with oscillatory (axisymmetric) solutions with two conjugate complex eigenvalues. In particular, the aforementioned 
mechanical examples can often be described by a Mathieu equation 



where x denotes an oscillating state, a is related to the fundamental frequency of the system and S is the amplitude 
of the periodic perturbation (e.g. pendulum with vertical oscillations of their pivot, fast spinning disks with time- 
varying loads [29[). However, the dynamo problem can be reduced to a Mathieu-like equation only in a few simplifying 
cases. For example, the thin disk approximation used in the galactic aw-model from [3lj allows a description of the 
amplitude of the magnetic vector potential by an equation similar to ((SJ) . A parametric resonance (also called swing 
excitation) is observed when the frequency of the perturbing velocity pattern is twice the natural frequency of the 
dynamo. Likewise, parametric resonances have been found in a Bullard-type disc dynamo model with a periodic 
modulation of the disc rotation. The model is based on a damped variant of {5} and resonances which allow magnetic 
field excitation for rather low velocities of the conducting disc are also observed at higher frequencies [48]. A more 
general description for systems subject to periodic distortions is given by 



where x is the departure from a basic state, M(t) is periodic with the forcing period T — 2it/uj and A is an operator 
that describes the unperturbed system. In principle, equation Q is suitable to describe the system examined in the 
present study (with A denoting the unperturbed dynamo operator). However, without detailed (analytical) knowledge 
of the eigenvalues and eigenvectors of the unperturbed problem an analytical treatment remains impossible, so that 
here we restrict ourselves to a qualitative description of the fundamental approach. The general form for the solution 
of system (|9|) follows from Floquet theory and is given by x(t) ~ J2i e Mi 'P j(t), where Pi(t) has the period T and 
Pi = e Mi * are the so called Floquet multipliers (e.g. j^Ho]). It can be shown that the computation of eigenvalues 
and eigenvectors of the problem yields different regimes that are characterized by stable, unstable, periodic or pseudo 
periodic solutions with simple and double periodicity in accordance with the frequency scaling in the modulated 
and/or resonant regime that is obtained in our study. Furthermore, a general calculation of the Floquet exponents \ii 
immediately explains the occurrence of exceptional points where eigenvalues and eigenvectors coincide. Concerning 
our present study, however, neither a reduction to a Mathieu-like equation nor an analytical computation of the 
Floquet exponents is possible without severe simplifications so that we abandon from further quantitative discussions 
at this stage. 

The observed behavior can also be explained on the basis of simple physical principles. For sufficiently slowly drifting 
vortices the system adjusts itself to an optimum state and the (azimuthal) phase between magnetic eigenmode and 
vortex pattern remains fixed so that the field growth becomes maximal. This state is essential characterized by 
the alignment between the magnetic eigenmode and the non-axisymmetric velocity mode which is roughly the same 
independently of symmetry breaking or vortex drift (see figure 1 1 3|) . Increasing uj v , the magnetic eigenmode cannot 
follow the ever faster (m = 2) vortex drift, but "bethinks" of its own eigenvalue (in the unperturbed state) to which 
it converges in the limit |w v | — > oo. In doing so, it will be "beaten" by the m = 2 vortex mode with an ever increasing 
frequency 2w v , which explains the occurrence of the second frequency involved. 

V. CONCLUSIONS 

We have examined kinematic dynamo action of a von-Karman like flow of a conducting fluid in a cylindrical 
container. When the flow breaks the ideal equatorial symmetry of the system, the critical magnetic Reynolds number 
for the onset of dynamo action increases with the amount of symmetry breaking and a time dependence is introduced 
in terms of an azimuthal drift motion of the dominant dynamo eigenmode. The frequency of this drift increases with 
the amount of symmetry breaking as well as with the magnetic Reynolds number. 

The main focus of our examinations has been on the interaction of this field drift with a non-axisymmetric time- 
dependent velocity perturbation and the resulting impact on dynamo action. In summary, what we observe is the 
following: the temporal behavior of the system is governed by three different time scales: the decay/growth time, 



x + (a — S cos(2t))x = 0, 



(8) 



x(t) = M(t)x(t) = (A + Bsmujt)x(t) 



(9) 
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FIG. 13. Alignment of magnetic eigenmode and velocity perturbation for different values of the vortex drift frequency cj v . 
All runs stem from the resonant regime. The color coded structure denotes B v and the contour lines show the axial velocity 
perturbation u v z . Left: a — 0, right: a = 0.62 



the magnetic field drift and the vortex drift which, in contrast to the first two represents an imposed quantity. For 
rather slow vortex drift frequencies, the first magnetic eigenmode (dominated by an (m = 1) component) is enslaved 
by the drifting vortex pattern, hence we see here Wf ~ uj v . Connected with the linear frequency relationship in this 
regime we observe a parabolic shape of the growth rate, with a maximum close to the resonance point where the 
vortex drift frequency would roughly correspond to the eigenfrequency of the magnetic field for the unperturbed, 
axisymmetric flow, and a quadratic reduction of the growth rate nearby this point. Outside of the resonant regime 
the field amplitude and the field drift are modulated with twice the frequency of the vortex drift. 

The simulations presented in this study show similarities with the results presented in [ToT ] where dynamo action 
was examined using a flow field field obtained from non-linear simulations of spherical s2t2 flow in a sphere. Kinematic 
dynamo simulations using the time averaged flow or different snapshots of the velocity field did not exhibit dynamo 
action, whereas this was indeed the case when considering the time-dependent flow. The exclusive occurrence of 
dynamo action with a time-dependent flow was interpreted as dynamo action based on non-normal growth originally 
described in (5lj . The essential prerequisite for non- normal growth is an appropriate mixing of (non-orthogonal) 
eigenstates through a time dependent linear operator which can lead to growing modes even if the contributing 
eigenstates alone correspond to decaying solutions in a stationary system. However, we observe a resonant behavior 
even in runs without equatorial symmetry breaking and with stationary vortices, i.e. in systems without any time 
dependence so that no time-dependent contribution is available that may provide for a mixture of non-normal modes. 
Hence, the role of non-normal growth in our models and the comparability with the model of [16| must for now remain 
an open question. 

It is difficult to conclude, if the resonance effect can be realized in existing dynamo experiments. Although a 
coincidence of forcing frequency (the frequency of the drift motion) and magnetic field frequency cannot be ruled 
out, the resonance condition is a quite particular case and most probably can only be fulfilled by chance. This is 
particularly true in the case of a more realistic non-linear analysis, in which the Navier-Stokes- and the induction 
equation are coupled so that the drift of the vortex pattern might strongly be influenced by back-reaction of the 
magnetic field. Furthermore, in the examined VKS-like configuration, the vortex drift frequencies observed in the 
water experiment are far away from the resonance condition (at least for reasonable values of the equatorial symmetry 
breaking) so that it is unlikely that the vortices affect the VKS dynamo in its actual set-up. Nevertheless, it might 
be suggestive to change the large scale flow geometry in order to adjust the relation of vortex drift and field drift, e.g. 
by changing the aspect ratio [52| or by fixing the vortices using some wall inhomogeneity. 

The occurrence of parametric resonances might also play a role for future laboratory d yna mos, like the large 
scale precession dynamo scheduled at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR) [53j ]. It is already well 
known from preliminary water experiments that drifting non-axisymmetric vortices aligned with the symmetry axis 



14 



of the cylindrical container emerge prior to an abrupt transition to turbulence 54J . In the appropriate regime, that 
means with suitable vortex drift frequencies and perturbation amplitudes, a parametric resonance similar to the effect 
examined in this study can be useful for increasing the efficiency of the dynamo process. Such effects indeed have 
been suggested for the Earth's core flow where, however, parametric instabilities caused by tidal forcing may provide 
the required periodic forcing [55| . 
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